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1 Introduction 

Let C M'^ be an open set. For any T < oo, set Qt = ^ x (0,T). The 
motion of a homogeneous, incompressible fluid through Q is governed by the 
following equations 

' dtU-divS +{V ■V)u + Vn = f, in Qt 
V-u = 0, in Qt (1-1) 

u\t=o = uo{x), in Q, 
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where u is the velocity, n is pressure and / is the force, V is chosen a solenoid 
vector function and tangential to the boundary of Q, Uq is initial velocity and 
S = {sij)^,j^i is stress tensor. The above system (1.1) has to be completed 
by boundary conditions except that Q is the whole space and by constitutive 
assumptions for the extra tensor. Concerning the former we can impose the 
following Navier slip boundary conditions 

u-n = 0,{S ■n)r - aur = 0, on dQ x (0, T), (1.2) 

where a is the frictional constant. 

Many extra tensors are characterized by Stoke's law 5* = iyD{u), where 
D{u) is the symmetric velocity gradient, i.e. 

Assume that u is a. constant and y = m, (1.1) is called incompressible Navier- 
Stokes equations. 

However, there are phenomena that can be described hj u = iy{\D{u)\) 
with power-law ansatz to model certain non- Newtonian behavior of the fluid 
flows, and they are frequently used engineering literature. We can refer the 
book by Bird, Armstrong and Hassager [15] and the survey paper due to 
Malek and Rajagopal [30]. Typical examples for this constitutive relations 
are 

SiDiu)) = f,i6+\D{u)\r'Diu) 
SiD{u))= fi{5+\D{u)\'')'^D{u), 

with 1 < g < oo, 5 > 0, and /i > 0. 

The mathematical analysis of these models started with the work of 
Ladyzhenska [3^, [3l], [35]. She investigated the well-posedness of the initial 
boundary value problem with non-slip boundary conditions, associated with 
the stress tensor (1.3). In 1969, J.L. Lions [SB] proved some existence results 
for Laplacian equation with p > 1 + -^^^ and the uniqueness for p > 
under no-slip boundary conditions. In those papers, the authors applied the 
properties of monotone operator and Minty trick theory for the stress tensor 
satisfies the strict monotonicity and coercivity. 

Over these years, Ladyzhenska's and Lions' work were improved in several 
directions by different authors. In particular, for the steady problem, there 
are several results proving existence of weak solution in bounded domain [22l 
[25l|26] , interior regularity [HE?] and very recently regularity up to boundary 
for the Dirichlet problem [5HTT | [T ^ [2CT | Hn]. Concerning the time-evolution 
Dirichlet problem in a 3D domain, J. Malek, J. Necas, and M. Ruzicka [2S] 
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study the weak solution for p > 2. Later, L. Diening et.al have recent 
advances on the existence of weak solutions in [22] for P > | and in [23] for 
j9 > |. There are also many papers dealing with regularity of for evolution 
Dirichlet boundary problems and we refer instance to [3tlH[8]-[TT | [T6 | [T7]. In 
the three-dimensional cube with space periodic boundary conditions, there 
are a lot of literatures for the well-posedness of this model, we refer to the 
monograph [2B] and papers [1^1^. 

It should be emphasized that theoretical contributions mostly concern the 
homogenous boundary condition and space periodic boundary conditions. 
However, many other boundary conditions are important for engineer exper- 
iment and computation science. Commonly used boundary conditions are 
Navier-type boundary conditions, which were introduced by Navier in [5S] . 
Newtonian fluid under Navier slip boundary conditions was studied by many 
mathematician, [T2|[T3] and [12]. However, there are not too many results 
for non- Newtonian fluid. In [5112^ , the authors investigated the regularity of 
steady flows with shear-dependent viscosity on the shp boundary conditions. 
M. Buh'cek, J.Malek and K.R. Rajagopal [IB] obtained the weak solution 
for the evolutionary generalized Navier-stokes-like system of pressure and 
shear-dependent viscosity on the Navier-type slip boundary conditions in 
the bounded domain. 

In this paper, we consider the problem (1.1) with stress tensor S induced 
by JO— potential as in Definition 2.1, when Q = M.^, under the following slip 
boundary conditions 

u ■ ^1,3=0 = 0, iiSiD{u)) • n) - (n ■ SiD{u)) ■ n)n) 1,3=0 = 0. (1.4) 

In fact, this problem corresponds to the free boundary problem for the non- 
Newtonian fluids with free surface supposed invariable. 

Since we choose the stress tensor induced by a p— potential, and then we 
will obtain the equivalent conditions: 

= 0(i=l,2). (1.5). 

From these conditions, we extend to the external force term / and initial 
velocity uo to whole space by mirror reflection method and change (1.1) 
into a Cauchy problem. Hence, we can focus on the regularity estimates, 
uniqueness and existence of this Cauchy problem. Then one can obtain the 
existence of the solution by Galerkin Method in the half space. 

The paper is organized as follows. In section 2, after recalling the no- 
tation and presenting some preliminary results, we give the definitions of 
the p— potential and weak solutions. We also present the existence of the 



Us 



dui 



X3=0 
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divergence-free base with boundary conditions (1.5) in W"^'"^. In section 3, 
we show some theorems for the existence, regularity and uniqueness of weak 
solutions for the system (1.1) with boundary conditions (1.5). 

2 Preliminaries 

In this section, we will give some assumptions, function spaces and definitions 
for weak solution. We will show the Korn-type inequalities for unbounded 
domain and construction of the basis with boundary conditions (1.6). Let 
j^nxn ^YiQ vector space of all symmetric n x n matrices ^ = We 
equip M"^" with scalar product C '■ V = Yl^j=i^ijVij norm |^| = (^ : 77)2. 

Definition 2.1. Let p > 1 and let F : M+ U{0} ^ U{0} be a convex 
function, which is on the ]R+U{0}; such that F(0) = 0, F'(0) = 0. 
Assume that the induced function $ : M"^" — > ]R_(_|J{0}, defined through 
$(S) = F(|S|), satisfies 

Y,i.d,kdim^){B)C,kCi^ > 7i(l + l^n^lCp, (2.1) 

jklm 

l(VL„$)(5)|<72(l + |Sr)"^' (2.2) 

for all B,C ^ M"^" with constants 71, 72 > 0. Such a function F , resp. 
is called a p— potential. 

We define the extra stress S induced by F, resp. $, by 

for all B e M"^" \ {0}. From (2.1), (2.2) and F'(0) = 0, it easy to know that 
S can be continuously extended by S{0) = 0. 

As in the |2T] and [2H], one can obtain from (2.1) and (2.2) the following 
properties of S. 

Theorem 2.2. There exist constants Ci,C2 > independent 0/71,72 such 
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that for all B,C e Af^^" there holds 

5(0) = 0, (2.3) 
5^(5,,(S) - S,,{C))iB,, - a,) > ci7i(l + \B\' + \C\'y-^\B - 

Y,S^,{B)B,, > ci7i(l + \B\^y-^\B\^, (2.4) 

\S{B) - S{C)\ < C272(l + \B\^ + \C\^)'^\B - C|, 

\S{B)\<C2-/2{l + \B\^)'^\B\. (2.5) 

In the following part of this section, we will give some function spaces 
and the definition of weak solutions for the system (1.1). 

D{Rl) = {ue C^{Ml) : = on = 0}, Vp{Rl) = {u G D{Rl) : V ■ u = 0}" 



Wlp 



p3 



H = Dijk^j 

Denote Qr = {x G R\ : < R} for R > , then we have the corresponding 
spaces for domain Qf> as follows, 

D{Qn) = {ue Co^(n^) : ^3 = on X3 = 0}, Vp{Qr) = {ti G D(l]«) : V ■ « = 0}"^ ""', 
Let T]^ = BrD {x3 = 0}, r| = OBr n Rl and Q = Rl x [0,T], = 

X [o,r]. 

Definition 2.3. Lei | < jo < oo, under the assumption of Definition 2.1. Let 
f E H or f E V* , which is the dual space of Vp, and uq E H with V ■ uq = 
in the sense of distribute. A vector function u G L°°(0, T; H) fl -L^(0, T; Vp) 
is called a weak solution to (1.1) if the following identity 



{u-dt<i))dxdt+ / {S{x,t,D{u)) -V ®u) : D{<p)dxdt 

(2.6) 

/ ■ ^dxdt + / uq- 0(O)dx 



holds for all (p E C°°(R% x [0,T]) with div(j) = 0, 031^3=0 = 0, and supp(p C 
X [0,T). 

We need to recall the following Korn-type inequality (see Theorem 3-2 
in [32].) 
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Theorem 2.4. Let K he cone in M" and p > 1. If Jj^ \ D(u)\'Pdx < +00, 
then there is a skew- symmetric matrix A with constant coefficients such that 

[ \V{u{x) - Ax)fdx <C [ \D{u)\Pdx 
Jk Jk 

where the constant C does not depend on u. 

The previous resuh leads to the following. 
Corollary 2.5. There exists a constant C depending only on p such that 



\Vu\\p < C\\D{u)\\p 



for all u E C( 



00 / 
I 



Proof. Since the domain is a special cone in M^, therefore, along the proof 
Corollary 1 in [27J, it is easy to get the result by Theorem 2.5. □ 

To construct the basis in with the boundary conditions (1.5), 

we consider the following problem 



' - Au + Vp = f, 
V ■u = 0, 

dui du2 



% = 0, 



dx-x dx-^ 



in VIr 
in VLr, 

on rjj, 



(2.7) 



^ M = on r^, 

The following definition 2.6, Lemma 2.7 and its proof will be found in 

Definition 2.6. By a weak solution of the problem (2.7) we mean a function 
u{x) G V2{^r) such that 

{D{u),D{v)) = {f,v), WveV2{nR) 

Lemma 2.7. Assume f G H{Qr), then there exists a unique solution {u{x) , p{x)) 
to problem (2.7) such that u G V2{^r) H W'^''^{VLr). Moreover, the following 
estimates hold: 



P(«)l|L^(n.)<C||/b(0,) 

||V^M||i2(f^^) + ||Vp||l2(j^^) < C + \\u\\v2{nn)) 

||V^/||L3(n^) < C 

where C is independent of u, f. 



1 j. 



(2.8) 
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With the aid of previous lemma, one can prove the following proposition 
Proposition 2.8. The eigenvalue problem 

— Au + Vp = \u, in 

V • M = 0, in Qr, 

dui du2 



uz = 0, 



w = 



on rl 



on r 



A e R, u & V2{flR) admits a denumberable positive eigenvalue {Aj} clustering 
at infinity. Moreover, the corresponding eigenfunctions {ai} are in W'^''^{Qr), 
and associate pressure fields Pi G W^'^{flR). Finally, {ai} are orthogonal and 
complete in H^Qr) and V2{flR). 

Proof. The mapping A : f — > u defined by Lemma 2.7 is linear and 
continuous from H{Qr) onto V2{^Ir), into W^'^{Qr). Since Qr is bounded, 
by Rellich Theorem, we know that W^'^{il,R) ^ L'^^Qr) is compact. It is easy 
to know that operator A is a positive symmetric and self-adjoint operator on 
L^ftR). Therefore, A possess an sequence of eigenfunctions af. 

Aui = Xitti k > 0, Xi > 0, Aj — > cxo as A; — )■ oo 
{ai,aj)L2 = 5ij, {D{ai),D{aj))L2 = XkSij. 

By Lemma 2.7, we can get for each i, there exists Pi with the estimates 
(2.8). □ 



3 Main results and their proofs 

To study the well-posedness of problem (1.1), we define a reflection as follows 

{{Ui{xi, X2, X3),U2{xi, X2, X3) , U^ixi, X2, X3)) if X3 > 0; 

{Ui{Xi,X2, -X3),U2{Xi,X2, -X3), -U3{Xi,X2, -X3), if X3 < 0. 

(3.1) 

Next, we will show the existence, uniqueness of strong solutions to the 
problem (1.1). We give the definition of the strong solution for the problem 
(1.1) as follows 
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Definition 3.1. We say a couple (m, tt) is a strong solution to problem (1-1) 
if 

u e L~(0, T; Wl^M) n i^^O, T; W^fM) ^ ^'(0, T- V,) n L-(0, T; //) 
^ e L2(0,T;LL(n));7r e L^'(0, T; Lf^^p. 

(3.2) 

where p' = and satisfies the weak formulation 

f ^ipdx + f S{D{u)) : D{ip)dx + f {V-V)u- cpdx 
Jq (Ji Jn Jn g-j 

= / Tvdivcpdx + / fcpdx. 
Jq Jq 

holds for all (p e C^(Q) anc? almost allt e (0, T), at same time, the boundary 
conditions hold in the sense of trace. 

At first, wc provide the definition of difference and recall a well-known 
result. Fixed any domain fl C M^,0' CC fl, and we put 6{fl',fl) — 
dist(l]',ai]\{x3 = 0}). 

Definition 3.2. For g : Q — > IR'^ we set 

{Ax,kg){x) = g{x + Xck) - g{x), xen',0<X< S{n', n),k = l---3. 

where €1,62, 63 is the canonical base ofM.^. We shall omit the dependence on 

k where the meaning is clear. 

Lemma 3.3. For any u e W^'P{n) and < |A| < 6{n',n) it holds 

\\Ax,ku\\p,n' < |A|||M,fc||p,f2- 

For above Q, Q' and 6{Q', Q), then the following theorem and lemma show 
that the regularity and uniqueness of the weak solutions to the problem (1.1). 

Tiieorem 3.4. Let | < p < 2, / G Lp'{Qt), V G L'^{0,T;W^'\n)), Uq G 
V2 nif satisfy the boundary conditions (1.5), and S be given by a p-potential 
from Definition 2.1. Ifu G i/(0, T; Vp) nL°°(0, T; H) is the weak solution for 
problem (1.1), then this solution is also a unique strong solution to problem 
(1.1) such that 

), (3.4) 

/ / (1 + \D{u)\Y-^\VD{u)\^dxdt < C{\n'\,uoJ,T, ||y|Uoo((o,T)xQ)). 

Jo JQ' 

(3.5) 

-g^\\l^(0,T)xQ') + \MD{u))\\l'^(0,T-,LHQ')) < C{S{Q' , , Uq, f ,T, \\V\\loo^^o,T)xQ)) ■ 



(3.6) 
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Proof. Firstly we extent the Uo and force term / to the whole space by 
the reflection defined in (3.1), however, since V G VF^'^, we need apply the 
extension of Theorem 5.19 in [2], Denote these functions by Uq, /*, V* 
respectively. 

We begin to consider the Cauchy problem as follows 

' dtv - diYS{D{v)) + {V* ■ V)v + Vtt = /*, in X (0, T) 

< V-v = 0, in X (0,T), (3.7) 

^ v\t=o = u*oix), in M^. 

From [39J, There exists a weak solution u G Lp(0, T; Vp{R^))nL°°{0, T; L'^{R^)) 
to problem (3.7), then by interpolation inequality, one follows that u G 
L^((0,r) X R^). If p > I then p' < f , and V* u e LP'{{0,t) x SJ, 
since V* G Along the proof of Theorem 2.6 in ^4jy, we know that 

lkllLp'((0,i)xBa) — ^(\\'^\\Lp'{{0,t)xBa) + 11*^ + /*llLP'((0,i)xBa) + ^) 

where Ba is any ball of R^ with radius a, C only depends on p, T, /, uq, a. 

For any p > such that < p < 5{il', il), Set fip = {x G fi; dist(x, fi') < 
p}. Now fix r < j6{Q',Q), there exists a ball Ba, such that C -Ba Let 
us choose a cut-off function r] such that = 1 on i^^, ''7 = in M'^ \ Q2r, 
< 7] < 1 and iVr^l < ^, |V^?7| < ^ in fi2r, where the constant C depends 
only on the geometry of dQ. 

If|A| < r, it results that A_A(??^AAn) G ^^((O, T) xfi3,)nLP(0, T; iyo'^((]3,)), 
but it is not divergence free. Take A_x{r)'^Axu) as a test function in the first 
equation of (3.7), we obtain 

{ut,A.x{v^Axu)) + {S{D{u)),D{A^x{v^Axu))) + {V* ■ Vu, A^xiv^Axu)) 
= {n,dw{A^x{v'^xu))) + {f,A^x{v'^xu)). 

Set 

Ji = {ut,A^xiv^Axu)), 
J2 = {S{D{u)),D{A^x{v'^xu))), 
J, = {V*-\/u,A_xiv'^xu)), 
Ji = {7i,dw{A^x{v^Axu))), 
,h = U\A^x{r,^Axu)). 

Clearly, Ji = f^\\r^Axu\\l,^^^^y Let Ix{u) = j^Jl + \D{u){x + Ae^)! + 
\D{u)\Y-^\r]AxD{u)\'^dx. 
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J2 = 2 /" S{D{u))A_xsym{Axu^r]Vr])dx- I rj\Ax{S{D{u)))AxD{u))dx 

J2I — J22- 



Since (2.3), 

J22>Ci7i I (l + |L>H(x + Aefe)| + |L>(M)|r->A;,L>(M)|Mx = Ci7i/AH, 

\J2i\ < c\\S{D{u))\\^p'^^^J\A_xsym{AxU(^r]Vr])\\LP(n^^), 
Thus 

\\A_xsym{Axu (g) r]Vr])\\LP(n3r) < \M\\^syM^xu(g) 7]Vr])\\LP{nsr), 

< |A| [\\Vr]\sym{Axu^Vr])\\\LP(n3r) + \\vsym{AxVu^^^^r})\\Lv{n3r) 
+ \\sym{Axu^Vr)^^)\\\Lv{nsr)] 

< 4C^||V^.|U.(03.) + ^ f / Iv'^iAxuWdx) " . 

^ ^ \Jn2r J 

Since rfQ{Axu) = V{r]Axu) - (Vr)) ■ Axu, thus 

/ \rjV(Axu)\Pdx) " < ^||Vti|U.(n3.) + cJ I \r^D{Axu)rd^ " 

JVl2r J r \jQ.2r ' 

However, by Holder's inequahty we have 

/ |?7D(AA«)rdxy <lx{u)^{\ {\-r\D{u){x-r \eu)\\\D{u)\fdx 
Hence 

J2 < \\S{D{umLp' insr) f ^l|V«|kp(n3.) + CIx{u)Hl + \\Vu\\%j) - C^Jx{u) 



(ri\2 ^ 2-p \ 

^||V«|Up(03.) + CIx{u)Hl + ||V«|| J(^3^)) j - CjJxiu) 

< \'Ci\nsr\,r,e) (1 + ||Vu|r^,(^3^)) - (Ci7i - e)/^^. 

\J,\ = \{f,A.x{v'^xu))\<\\f\\ 

LP{n2r) 

< C{r,p)X^f\\^p>^^J\Vuhpin,.) + \X\\\v'AxVuhpin2.) 

< \X\Ix{u)H^ + l|Vu||5j'n3.))ll/llL.'(^3.) + ^(^ + l)^"\\f\\Lp'insr)WUns.) 

< Cip, \Qsr\)X' (ll/ir;.'(^3^) + l|Vu||I,(^3^) + 1) + elxiu). 
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Prom the estimate for pressure, divergence-free and the method above, we 
have 

\Ja\<2 



TTAx{r]r]iAxUi)dx 
d 

< C I A 1 1 1 TT 1 1 (Ha.) \\-Q^{VVxAxUi)\\Lv (nsr) 

1 CI Al 

< (lkr;.'(n3.) + l|V-lli.(n3.) + l) + ^h{u). 
Now we estimate the term J3. In fact, 

J3 = / v'^AxV* ■ VuAxudx -2 riV* ■ Vr]\Axu\^dx 

'■— J3I + <^32 



Since | < p < 2, then 2 < q = < ^ = p*, by Holder inequahty, we 
have 



J,^ = \\r]AxV*\\ 

<C|A|||Vti|U.(n3.)l|VV*||La(f73.)II^A,ti|U,(n30 

< C|A|||Vn||LP(f23,)||l^*|U2,2(f^3^)||r/AAti||L«(Q3,) 

< C'|A|||y||ic^((o,T)xn)||Vu||LP(n3.)||^AAM||L«(n3,) 

C 

J32 < -\\V*\\Le(n,^)\\r}Axu\\Li(n3r)\\^xu\\LP{nsr) 

< C|A|||1/||L«=((o,r)xQ)||Vu||LP(n3.)||r?AAw||L,(Q3,). 



It is imphes 



iJsl < C(J, \n3r\)\M\\V\\L^{{o,T)xn)\\'^u\\LP{nsJ\v^xu\\L'>{nsr) 
Since 2 < q < p*, from the following interpolation inequahties 

IUIIl'j ^ lUlliP* IUIIlz^; 

where 9^^,01^ We obtain 
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where 

7 12 

Qi = (1- 6)a = — a, < a < 1, 

5p — 6 

Q^ = ae + {1- aOi) = + ^^^(1 - a) 

op — 6 2p 

7n — 12 

Q3 = (l-«)(l-ei) = ^^^(l-a). 

Since 

II^A,«llL.*(n..) < C||V(r/A,«)|U.(n,.) < C||L»(r/A,«)|U.(n,.) < Ch{u)Hl+\Vu\%^^^). 
It infers that 

By Young's Inequality, wc get 

\Js\ < C^||V^IU^((o,T)xf.)|A|(^+«^)^' (l + llVx.llit^S'^'''^'') lkA,t.|e4^)+e/,H. 
Then choose 5 and 5' satisfy the following identities 

^S^l, (l + Q, + l^Q,)5'^p, and^ + i = l. 

From these identities, we can obtain a = -^^^^^i^-^. Since p > |, thus 

< a < 1, (1 + Qi + Q3)6' = 2 and ^ < 1. Therefore, 

1 Jsl < C||y|Uoc((o,T)xf.)|A|(^+«^)^' (l + ||Vx.||i,(^3^)) |kA,x.|g4^) + eh{u). 
Combined these relations of Ji, • • • , J5, we can conclude that 

^^Wv^xuWh^^^^) + (Ci7i - ^e)h{u) 

< cx^ (1 + iiv.r,.(,3^) + wmC^,^^^ + wfCJ (3.8) 

+ C||y|Uoo((o,T)xf.)|A|(^+«3)^' (1 + ||V^.||i,(^3^)) ||r/A,u||?4^) 
Since (1 + Qi + Qs)^' = 2 and set 

r{t) = C||V-|Uc.((o,T)xf2) (1 + l|Vu||i.(^3^)) , 

(3 9) 

Mi) = c (1 + iiv.r,.(,3^) + ikr;,^,^^^ + mCJ^ 
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where C does not depend on A, u. Hence (3.8) can be rewritten 

^ll^llL^(n..) + (^^i7i-46)^ 

<Mt) + r(t)||^l|?4^, 

Since p > |,it is easy to know that < 1. Prom Young's inequahty, we 
have 

^ll^llL^(n..) + (^?i7i-46)^ 

<(Mt)+r(t))+r(0||^||i.(^^^). 
By Gronwall's inequahty, it foUows 



< 



\\h{n,r)(^) + h{s)exp(^- r(r)dr^ds exp^ r(s)ds. 



A 

Assume that Uq G l^, then it imphes that for aU t G [0,T] 



I^IIW,W + (C,7.-4e)/Mf) 



< 



exp / r{s)ds. 



(3.10) 

Choose e = ^ and from (3.9), (3.10), we conclude for any t & [0, T] 

l|V«||i.(n.)(t) + fl{u)<C. (3.11) 

where/(w) = /^^(l+2|D(M)|)^'-2|i:)(VM)pda; and C depends on \\V\\Loo^(^o,T)xQ),P,a,Uo, f,T,r. 
Since 

for aU i e [0,T], thus ||w||LP(o,r;VF2.p(f2,.)) ^ C*- Multiply the first equation of 
(3.7) by rj'^Ut and integrate on the $13^, one obtains 

h«t||i^(n3.)(^) + ^ll^'*(^(«))llLi(n3.) 

<2 I \VriS{D{u))riUt\dx+ I V* ■Vurj'^Utdx (3.12) 
+ 2 / ttVt/ • rjutdx + f • rfutdx. 
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Since ^ + 4 < 1, by Holder inequality, one can 



p 



\Wr]S{D{u))r]Ut\dx <\\S/r]\\L, (^3,.) \\S{D{u))\\^^> 1 1 I in^r) 



/ V* ■ Vurfutdx < C||y||Loo((0,T)xn)hVli||L2(f^2r)ll^^jL2(f^3,,) 

Thanks to Young's inequality, we can obtain that 

hWt|lL2(O3,x(0,t)) + \W^{D{u))\\L^(o,t;L^{n;,)) < K 

where K depends only on p,a,uo,f,r, \\V\\Lo^^(o^T)xn), || Vw||lp(r3x(o,t))- 

In fact, from this proof, we can see that the bound depends on the 
measure of Q and Q'. Hence, if the radius of the ball B is fixed, then 
II V^m||lp(o,t;Lpb) < C, consequently, u G C^i^B). We use the following argu- 
ment ( see [27]) to know that u{x, t) — )■ for almostt G (0, T), as |a;| — )■ 00. Let 
the radius B be one, suppose that there exist e > and a sequence {x„} C M'^ 
with lim„^oo l^^nl 00, such that for almost t G (0,T), u{xn,t) > e. By the 
continuity of u{x,t), then we get that if \x — x„| < S = min{l, {^)~} then 
u{x,t) > |. Without loss of generality, we assume that \xi — Xj\ > 2 provided 
that i ^ j thus 

/ ( / \ufdx)'^dt [ if |Mr*dx)^dt = +00. 

Jo j Jo Jb{xj) 

and this contradict with the fact u G L^{0,T; ^^(IR^)). 

From the fact u{x,t) — )■ a.e. t G (0,T), as |x| — )■ 00 and estimate (3.13), 
we can conclude that the weak solution is unique. Actually, assume that u, v 
are that weak solutions of problem (3.8), then set e = u — f and multiply the 
difference of the equations of u and v by e. After integrate over M^, we have 

{dte, e) + {S{D{u)) - S{D{v)),D{u) - D{v)) = 0. 

It reduces to 

^I|e||i2(i,3) + {S{D{u)) - S{D{v)), D{u) - D{v)) = (3.13) 
Use Gronwall's inequality in (3.14), we have e = 0, i.e. u = v. Define 

U* = {Ui{Xi, X2, -Xg), U2{Xi, X2, -Xg), -^(Xi, X2, -X3) 

71* = 7r(xi, X2, —X3), Vx G M^. 
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Then by the method in [13], the couple {u*, n*) is also a solution to problem 
(3.4) a.e. in M.^ x (0,T). Hence by the uniqueness, we know that u* = u. 

From the regularity of u, and Theorem 7.1 of [13] and Sobolev imbedding 
theorems, we know the solution u{t) to problem (1.1) is simply the restriction 
of u*{t) to the half-space and || VM*||Lp(]K3x(o,r)) < C\\^u\\lp(k.\x{o,t)), 
||'W*||l°°(o,T;L2(r3)) < C||'u||loo(o,t;H), and u{x) is a continuous function on 
Br for any r > 0. So that from Uz{x) = —U3{xi,X2,—X3) we know that 
11313.3=0 = 0. Analogously Ui{x) = Ui{xi,X2, —x^), Wx G M.^{i = 1,2), satisfy 
the conditions (1.5) in the sense of trace. 

The theorem is completely proved. □ 

By the minor modification of the proof in theorem above, we can obtain 
the regularity results in the case p >2. 

Lemma 3.5. Let p > 2, f e W' {Qt) n L'^{Qt), uq e V2 n H satisfies the 
boundary conditions (1-6) , and S given by a p-potential from Definition 1. 
If u & L''^{0,T;Vp) n L°°{0,T; H) is the weak solution for problem (1.1), then 
this solution is also a unique strong solution to problem (1.1) such that 

ll'w||L°°(0,r;H/l.2(n'))nL2(0,T;iy2,2(n')) < C{5{yi\ Vl), l^^'l, Mo, /, T, ||V^||Loo((o_T)xn)), 

f [ (l + |D(M)|)f-2|VD(M)|Ma;dt <C(5(fi',fi),|fi'Uo,/,T,||y|Uo.((o,T)xn)). 
Jo Jw 




L2((0,T)xf7') + ll'^(^(w))l|L-(0,T;Li(n')) < C{S{Q',n), \n'\,Uo,f,T, || F|| ((o,T) xH) ) • 



lk||L2((0,r)xn') < C{6{il',il), \il'\,Uo, f,T, ||F||L°°((0,T)xf7))- 

Remark 3.6. Since divV = 0, thus if u is a weak solution of problem (1-1), 
then we can obtain the following estimates 

II'"IIl°°(0,T;H) + \\'^\\lp{0,T;Vp) - ^HI^oWh + II / II (R^) X (0,T) ^ " i^'^^) 

It is easy to see that C does not depend on a,r, \Q'\,uo,f, || V"||ioo(-(o,T)xf2) ■ 

From these regularity estimates presented in Theorem 3.4 and lemma 3.5, 
we can obtain the existence of unique weak solution to problem (1.1) stated 
by the following theorem. 

Theorem 3.7. Letp > §,/ G Lp' {QT)nL^{QT),V e L°^{{0,T),W^^^{Rl)),uo G 
V n H, and S given by a p-potential from Definition 2.1. Then there exists 
a unique weak solution for problem (1.1), u & L^{0, T; Vp) fl L°°(0, T; H) and 
satisfies the inequality (3.14)- 
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Proof. From the proof in Theorem 3.4, it is easy to see that the weak solution 
is unique. We will use standard Galerkin method to prove its existence. 
Let 



\JQr = \J{xeRl:\x\<R}. 



R=l 



R=l 



Fix R > 0, we consider the auxiliary problem (2.7) for the initial Uq = 
P{Xnii{x)uQ{x)) and external force term = xuni^)!- 

As Lemma 3.5 in [31], we construct a sequence of solenoidal vector func- 
tions Vr defined in Qr such that Vr ^ V in Lp(0,T; as i? oo, 
||^||L°°(o,r;vy2,2(j^^)) < C, here C depends only on the Geometry of Qr, but 
does not depend on the measure of Qr. We also know that — )■ -uq in ^2 H if 
and/^ / in Lp'{Q), as i? oo. 

Choose the sequence {a^} is the eigenvector of the operator A as in 
Proposition 2.8, then {af } is a basis W'^''^{^Ir) fl V2{^Ir). We look for the 
weak solution to the following problem 

' dtu - diYS{D{u)) + {u ■ V)u + Vtt = /, in Qr x (0, T) 



V - M = 0, 



Us = 0, 



dui du2 



dxR dxs 



in ^R X (0,T), 
on r^x (0,T), 



u = 
^ w|t=o = uo{x), 

We find the approximation solutions with the form 



on r^x (0,T), 
in VLr. 



k=l 



For simplicity, in the clear meaning setting, we omit the superscript R. 
Therefore, Ck^mif) solve the following system of ordinary differential equa- 
tions 



— ttfe) + {S{D{um)),D{ak)) - (V ^ Um, Vak) = (/, a^). (3.15) 
at 

Due to the continuity of S, V, the local-in-time existence follows from Caratheodory 
theory. The global-in-time existence will be established by the following a- 
priori estimates. 
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Multiply the equations (3.15) by Ck^m, then sum over k and integrate on 
(0,t). We easily obtain 



I |i2 



^+ / {S{D{Um)),D{Um)) = {f,Um) + \\uo\\mnn)- (3-16) 
^0 



Hence 

SUp^\\Um\\l2^^^){t) + \\Um\\lp^o,T;VAnn)) - ll^oUlf + ||/||^p'(K3)x(0,r) 

(3.17) 

FYom (3.15) and (3.17), we infer that ||M^nllLp'(o,r;(yp(Q«))*) < C{R,M). 
Prom these estimates, and Aubin-Lions Lemma, we have 

u'^^u', weakly in L^''(0, T; Fp(fi^)*); (3.18) 

Um^u, weakly* in L°°(0, T; /.^(fii?)); (3.19) 

weakly in (0, T; ^^(fi^)), ; (3.20) 

«^ ^ u, strongly in L'^(0, T; L'^(Oi^)) g G [1, ^); (3.21) 

S{D{u^)) S^, weakly in L^' (Qr x (0,r)). (3.22) 

It is easy to see that 

( ut, au ) + (^^, D{ak)) - {V ■ Va,, u) = (/, a^). (3.23) 

Multiplying both sides of (3.23) by Ck,m and summing over k we find 

{ut,Um) + {S'',D{Um)) - {V ■ VUm,u) = (/,«„). (3.24) 

Let us pass to the limit for m — ?> oo in to (3.24). By the convergence prop- 
erties (3.18), (3.20) and (3.22), we know that as m ^ oo 

[ {ut,Um) ^ [ = ||M||i2(Q^) - ||Mo||i2(n^), 

JO Jo 

Jo Jo Jo Jo 



Since 



(V ■ VUm, u) - (V ■ Vu, u) = {Ur,(e) ® u) ■ V(m^ - u)dx 



From (3.20) and (3.21), we know that V ®u E W' {^r x (0,r)) whenever 
p > |. Hence {V ■ Vumju) — > as m — )■ oo, since {V ■ 'Vu,u) = 0. 
Subtracting (3.24) by (3.15) and passing to limit as m — )■ oo, we get 

hm r{S{D{u^)),D{um))dt= r{S^,D{u))dt (3.25) 
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By the monotonicity property (2.4), we can write the following inequality 
f{S{D{ul)) - S{D{^)), (Diurr^) - Dimdt > 0, (3.26) 

For any G L^(0,T; Vp{Q)). Thus, pass to the limit as m goes to infinity 
into this relation and using (3.20), (3.22)and (3.25), we have 

/ {S^- S{D{^)),D{u)- D{^))dt>0. (2.27) 
Jo 

For aU ^ e Lp(0, T; Vp{n)), Take ^ = u - e^?, e > and ^ G Lp(0, T; Vp{n)), 
we have that ^ 

/ {S^ - S{D{u-eip)),D{(p))dt>Q, 
Jo 

Letting e — > and using the continuity of 5", we arrive at 

/ {S^ - S{D{u)), D{ip))dt > 0, e (0, T; Vp{n)) 
Jo 

Choose —(f in place of ip, we get 

/ {S^ - 5(D(u)), D(^))dt < 0, e LP{0, T- Vp{n)) 
Jo 

This implies that — S{D{u)) a.e. VLr x (0, T). Thus the existence of weak 
solution to problem above is proved. 

Next we must consider the limits as R tend to oo. Now we choose a 
sequence of real number {i?jv : G N} increasing to infinity. We set mat = 
^-Riv extend Mat to zero outside VLr^^ to obtain a function still denote 
un & L°°(0, T; H) n L^'(0, T; V^) and satisfies the following apriori estimate 

^SUP^ \WN\\H{t) + ||Miv||LP(0,T;y^) < hofn + ll/llip'(R3)x(0,T) (3.28) 

Clearly, M is independent of N. From (3.28) and interpolation inequality 

we have ||V/v ® '^n\\ ^v' [m? x(ot)) 

< C(M). let Y = Vp, Therefore, from the 

equations (3.23) we obtain 

\W'n\\lp'{0,T;Y*) < \\S{D{un))\\lp'(Q) + \\Vn ^ Un\\l^Q) < C{M). 

By the Aubin-Lions Lemma, we know that there exists a subsequence un^ ^ 
u in i/(Qi? X (0, T)), thus un^ u a.e. in x [0, T]. 
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We choose e Co°(^+ ^ [0^ ^)) with div0 = 0, 03|x3=o = 0, and supp0 C 
R\ X [0, T). There exists a number K = K{(f)) > such that G = supp0 ^ 
Q.Nk ^ [Oj ^) for k > K. Prom the formula (2.6), we can obtain the following 
identity 



(3.29) 



■dt(j))dxdt+ / S{x,t, D{un^)) : D{(j))dxdt 
JQ 

I (VATfc ® matJ : D{(t))dxdt = / • 0(O)da;. 

jQ Jffii 



From the estimate (3.28), we know that mat^. — ^ in L^((5), S'(x, t, D(M7Vfe) - 
S in the space L^' {Q) and jV^r^ (g) itArjL2(Q) < C'||wo||i?||^||L°°(G)- By Vitah's 
theorem, we have as A; ^ oo 

/ (matj, ■ dt(f))dxdt {u • dt(j))dxdt, 
JQ JQ 

/ (VWfc ® MivJ : D{(t))dxdt / (Viv, ® «) : D(0)da:dt. 
7q Jq 

Prom these convergence and the formula (3.29), we know that 

{S{x,t,D{uN,)) : D{(j))dxdt [ (S : D{(f))dxdt. 



Q 



One must check that S{x,t,D{u)) = S a.e. in R\ x [0,r]. It suffices to 
prove that 



{S{x,t,D{uNj) : D{(l))dxdt / {S{D{u)) : D{(j))dxdt. 
Q Jq 

Indeed, observing that for any > Rnk the solutions un satisfies the hy- 
potheses of Theorem 3.4 and lemma 3.5 with Q — Q^j^ and fixed a domain Q' 
such that GcQ'x (0, T) ccQx (0, T), we get that un e i/(0, T; VV^'PiQ')) 
and 



k7vJU2((0,T)xn') + ||VMArJ|Loo(o,T;L2(n')) 

\.iin. f.T.\\V\\Toarr<\). in the case ^ 



9 

+ ||V^M7vJ|i,p(Q/x(o,T)) < (^(^(fi',!]), |l]'|,Mo,/,r, ||V^||loo(g)), in the case -<P 



\\u'Nj\L^{(0,T)xn') + ||VMArJ|ioo(o,T;L2(n')) 

+ \\V'^UNj\L^n'x{o,T)) <C{6{n',n),\n'\,uo,f,T,\\V\\Loo(^G)), in the case p > 2. 
Prom the boundedness above and the Aubin-Lions lemma we obtain that 
Vun^ in LP{n' X (0,T)), Vun^ a.e. in G. 
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Therefore, S{x,t, D{uni^)) — >■ S{x,t, D{u)) a.e. in G, then by Vitah's theo- 
rem, we obtain as A; — >■ oo 

/ {{S{x,t,D{u^^)) - {S{D{u))) : D{<P)dxdt = 
jQ 

f {{S{x,t,D{unJ) - {S{D{u))) : D{cp)dxdt ^ 0. 
Jg 

Whence, this proves the theorem. □ 
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